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Resolved 7^ in hard collisions of virtual photons: QCD effects 



o 
o 
o 

(N 

-t— > 
o 

O 



Jiff Chyla and Marek Tasevsky a 

Institute of Physics of the Academy of Sciences of the Czech Republic 
Na Slovance 2, Prague 8, Czech Republic 

Received: date / Revised version: date 

Abstract. The manifestations of QCD effects on quark and gluon distribution functions of longitudinally 
polarized virtual photons involved in hard collisions are investigated. It is shown that for moderate photon 
virtualities and in the kinematical region accessible at HERA and LEP these effects are sizable and sig- 
nificantly enhance theoretical predictions based on contributions of transversally polarized virtual photon 
only. 
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1 Introduction 

In QED quantized in covariant gauge, longitudinally po- 
larized on-shell photons are present, but due to gauge in- 
variance decouple, order by order in perturbation theory, 
in expressions for physical quantities. For the virtual pho- 
ton with virtuality F] P 2 its longitudinal polarization, de- 
noted 72, does contribute to physical quantities and gauge 
invariance merely requires that these contributions vanish 
as P 2 — > 0. In a previous publication [0 we have dis- 
cussed the contributions of 72 to two physical quantities 
using purely QED formula for quark distribution functions 
of 7£ . In this paper we continue our investigation of the 
relevance of 72 in hard collisions by incorporating the ef- 
fects of QCD radiation on parton distribution functions 
(PDF) of 72 recently derived in In the next Section 
the rationale for introducing the concept of the structure 
of virtual photon is recalled, followed in Sections 3 and 
4 by a short review of the QED and QCD formulae for 
corresponding PDF. The numerical relevance of the con- 
tributions of resolved 72 with QCD improved PDF are 
discussed in Sections 5 for the LO and in Section 6 for the 
NLO QCD calculations. 



2 Virtual photon and its "structure" 

Let us briefly recall the virtue of extending the concept of 
partonic "structure" to virtual photons |||| : 

— In principle, the concept of partonic structure of vir- 
tual photons can be dispensed with as higher order 



a Work done within the Center for Particle Physics under 
the project LN00A006 of the Ministry of Education of the 
Czech Republic. 

1 In this paper the virtuality of a particle with four- 
momentum k and mass m is defined as I k 2 — m 2 I . 



QCD corrections to cross sections of processes involv- 
ing virtual photons in the initial state are well-defined 
and finite even for massless partons. 
— In practice, however, the concept of resolved virtual 
photon is extraordinarily useful as it allows us to in- 
clude the resummation of higher order QCD effects 
that come from physically well-understood region of 
(almost) parallel emission of partons off the quark or 
antiquark coming from the primary 7* — > qq splitting. 

For the virtual photon, as opposed to the real one, its 
PDF ^| can therefore be regarded as "merely" describing 
higher order perturbative effects and not its "genuine" 
structure. Although this distinction between the content 
of PDF of real and virtual photons exists, it does not affect 
the extraordinary phenomenological usefulness of PDF of 
the virtual photon. As shown in || the nontrivial part 
of the contributions of resolved transverse virtual photon 
(7^.) to NLO calculations of dijet production at HERA 
is large and affects significantly the conclusions of phe- 
nomenological analyses of existing experimental data. 



3 PDF of 72 in QED 

Most of the present knowledge of the structure of the pho- 
ton comes from experiments at ep and e + e _ colliders, 
where the incoming leptons act as sources of transverse 
and longitudinal virtual photons of virtuality P 2 and mo- 
mentum fraction y. To order a their respective uninte- 
grated fluxes are given as 



a /l + (l-y) 2 ) 1 2m 2 c y 
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More precisely their pointlike parts. 
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Phenomenological analyses of interactions of virtual pho- 
tons and their PDF have so far concentrated on its trans- 
verse polarization. Neglecting longitudinal photons is a 
good approximation for y — ► 1, where / 7i (y, P 2 ) — > 0, as 
well as for small virtualities P 2 , where PDF of 7^ vanish 
by gauge invariance. But how small is "small" in fact? For 
instance, should we take into account the contribution of 
72 to jet cross-section in the region Et > 5 GeV, P 2 > 1 
GeV 2 , where most of the data on virtual photons obtained 
in ep collisions at HERA come from? The present paper 
is devoted to addressing this and related questions. 

In pure QED and to order a the probability of finding 
inside 7^ of virtuality P 2 a quark with mass m q , charge 
e q , momentum fraction x and virtuality r < M 2 , is given, 
in units of 3e 2 a/2n, as M 



q<}™(x,mlP 2 ,M 2 ) 



4x 2 (l-x)P 2 



1 - 



M 2 



(3) 



where r lnin = xP 2 + m 2 /(l — x). The quantity defined 
in (||) has a clear physical interpretation: it describes the 
flux of quarks that are almost collinear with the incoming 
photon and "live" longer than 1/M. For r miu < M 2 the 
expression (0) simplifies to 



q^ ED (x,m 2 ,P 2 ,M 2 ) 



4x 2 (l-x)P 2 



xP 2 + m 2 /(l — x) ' 
which for x(l — x)P 2 m 2 further reduces to 

QED / 



q^ u {x,0,P\M*)=Ax{l-x). 



(4) 



whereas for x(l — x)P 2 «ra^ 
q ? EB (x,m 2 ,P 2 ,M 2 )- 



T 4x 2 {l -xf 



demonstrating the fact that in QED the onset of 7^ is 
governed by the ratio P 2 /m 2 . 



4 QCD improved PDF of 7^ 

QCD improved PDF of 7^ have been derived in the leading- 
logarithmic approximation and for 1 < P 2 <C M 2 in 
0. By "leading-log" we mean resummation of the terms 
(cv s lnM 2 ) fe at each order k of perturbative QCD. Note 
that for 7^, there is one power of In M 2 more at each or- 
der of a s , the additional one coming from the primary 
QED 7* — > qq splitting. In the case of 7^ the analogous 
splitting gives rise to the term (jij), which is constant in 
P 2 . The resulting expressions ^| exhibit typical hadronic 
form of scale dependence and contain /Iqcd as the only 
free parameter. The condition P 2 <C M 2 guarantees clear 
physical meaning of the resulting quark and gluon dis- 
tribution functions. Moreover, by staying away from the 



3 The parameterization of PDF of 7^ can be obtained from 
chyla@fzu.cz. 



region P 2 ~ M 2 we avoid the region where power correc- 
tions of the type P 2 I'M 2 are essential and, in fact, more 
important than the effects described by PDF. The restric- 
tion from below 1 GeV 2 < P 2 ensures that hadronic parts 
of PDF of 72 , which have not been taken into account in 
the derivation in Q , can be safely neglected with respect 
to the pointlike ones 

The relevance of resolved 7^ in hard collisions of vir- 
tual photons f\ depends on the theoretical framework one 
works in. In the next two Sections we shall discuss the 
effects of including resolved 7^ within the LO as well as 
NLO QCD calculations. The difference between the nu- 
merical relevance of resolved 72 in these two cases arises 
from the fact that parton level calculations contain at the 
order ev 2 a 2 some of the effects that go into the definitions 
of quark distribution function of 7J and n fj j . 



5 Resolved 72 in LO QCD calculations 

5.1 DIS on 7* 

In LO QCD the structure function K 7 of the virtual pho- 
ton is given in terms of quark distribution functions by 
the same expression as for hadrons ^ 



F](x,P 2 ,Q 2 ) 



E 



2xej (q l (x,P 2 ,Q 2 ) + q l (x,P 2 ,Q 2 )) 



4 This claim is based on experience with SaS sets of parame- 
terizations Q and the assumption that hadronic parts of PDF 
of 7£ can be related to those of 7J [Q] . 

7x, contributes to soft collisions and related quantities, like 
otot(7*p), as well, but we restrict our discussion to hard colli- 
sions. For the former, the reader is referred to ||. 

6 We disregard the consequences of the reformulation of QCD 
analysis of F£ proposed by one of us in j?| as they do not 
concern the main point of our discussion. 
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1. Comparison of the contributions of 7J and 7^ to Fj g 
and D eff (right) for P 2 =2,5 GeV 2 and M 2 = 100 GeV 2 . 
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In all existing phenomenological analyses only target "f T 
has been taken into account, despite the fact that for 
P 2 <C Q 2 experiments at LEP actually measure [] 

the "effective" structure function 

F^(x, P 2 ,Q 2 ) = F^x, P 2 ,Q 2 ) + ^ L (x, P 2 , Q 2 ) 

given as the sum of contributions from target 7^ and t£. 
In Fig. |a,c we compare, for two pairs of P 2 and Q 2 typical 
for current experiments at LEP, i 7 ^ obtained with SaSID 
parameterization Q of PDF of j T with the contributions 
from target 7^ evaluated using both the QED and QCD 
expressions for Ql{x, P 2 , M 2 ) discussed in the preceding 
two Sections. The contributions from <?£ D peak around 
x ~ 0.7, with QCD effects suppressing them at large x and 
enhancing them on the other hand for x < 0.4. The pres- 
ence of the term proportional to In M 2 in the expression 
for qx in both QED and QCD implies the dominance of 
j T at large M 2 , but one would have to go to very large M 2 
for 7£ to become negligible with respect to j T . For fixed 
M 2 the relative importance of 7^ with respect to j T grows 
with P 2 , but to retain clear physical meaning of PDF we 
stay throughout this paper in the region P 2 <C M 2 . 



5.2 Dijet production in ep collisions 

The measurement of dijet production in ep collisions pro- 
vides another way of investigating interactions of virtual 
photons [jl0],[ll]]. In general the corresponding cross sec- 
tions are given as sums of contributions of all possible 
parton level subprocess. The simplest way of demonstrat- 
ing the importance of contributions of resolved 7^ employs 
the approximation [jl2| in which dijet cross sections are 
expressed in terms of a single effective parton distribution 
function of the photon (either 7^ or 7^) defined as 

D cB {x,P 2 ,M 2 ) EE 

M*, P\ M 2 ) + 5.(3, P 2 ,M 2 )) + -G(x, P 2 ,M 2 ), 

i=l 

where the factorization scale M is conventionally identi- 
fied with (a multiple of) jet E T : M = kE t . In Figs. |]b,d 
the contributions to D c g from j T and 72 are compared for 
two pairs of P 2 and M 2 typical for HERA experiments. 
In addition to effects at large x, which are similar to those 
for F^ s , D e s gets a sizable contribution from 7^ at small 
x, coming from its gluon content. The rise of D c g at small 
x is particularly clear effect of QCD improved PDF of 7^. 
After this estimate, we now proceed to discuss the contri- 
butions of 7^ to dijet cross sections evaluated with HER- 
WIG 5.9 event generator at the parton level. We could 
have used for this purpose also JETVIP [jl3|, which we 
shall use later at the NLO, but using HERWIG at the LO 
allows us to 

— estimate hadronization effects, 

7 Neglecting the difference of the fluxes (|l|j2|), which is a good 
approximation at small y, typical for LEP experiments. 



— cross-check the modifications implemented in JETVIP 
in order to include the effects of 7^ . 

For the purpose of this study we have modified standard 
HERWIG 5.9 by adding the option of generating the flux 
of 72 combined with the call to QED or QCD improved 
PDF of 72. For j T the SaSID PDF were used. All calcula- 
tions were performed for 0.05 < y < 0.95, three windows 
of P 2 : 1.4 < P 2 < 2.4 GeV 2 , 2.4 < P 2 < 4.4 GeV 2 and 
4.4 < P 2 < 1 GeV 2 and the following cuts on parton E T 

E { t\e {2) > E T , E T = 5, 10 GeV. 

The effects of HI and ZEUS detector acceptances have 
been approximately taken into account by performing all 
calculations without any restriction on parton pseudora- 
pidity as well as for — 3 < 77 < 0. 

The results for the first window in P 2 and without the 
cuts on rj are presented as functions of 77, x 7 and Et in 
Fig. |[ The characteristic dependence of the contributions 
of resolved 7^ on y is illustrated by plotting for each of 
the distributions in 77, Et and x 1 also its ratio to that 
of 7y for the whole interval 0.05 < y < 0.95, as well as 
for three indicated subintervals. Except for x 1 close to 1, 
QCD improved PDF of j1 enhance its contributions to di- 
jet cross sections compared to those based on the purely 
QED. For y < 0.5 and x 7 close to or 77 ~ 2.5, the con- 
tributions of resolved 7^ amount to about 80% of those of 
7y, whereas on average this number is around 50%. Re- 
ducing the range of 77 to —3 < 77 < affects (see Fig. ||) 
mainly the distribution dcr/dx 7 by suppressing it at both 
endpoints x 1 = and x 1 = 1. The ratia of the contribu- 
tions of 72 and "f T are, however, affected only little by this 
cut. Increasing the photon virtuality enhances, as shown 
in Fig. |J, the relative importance of resolved 72 with re- 
spect to j T On the contrary, rising the threshold E%, from 
5 GeV to f GeV reduces it, as illustrated in Fig. ||, by a 
factor of about 2, since large Et require large x 7 , where 
quarks from 7J dominate. 

Summarizing the message of Figs. ^-|[ we conclude 
that in the region yl 2 <C P 2 M 2 w E\ the contributions 
of 72 are substantial, particularly for 

— small y, 

— low Et-, 

— x 7 < 0.5, corresponding to rj close to the upper edge. 

The cuts enforced by HI and ZEUS acceptances reduce 
the sensitivity to 72, but its contributions still make up 
typically 30 — 50% of those of 7^ and can be identified by 
their characteristic y and P 2 dependencies. 



6 Resolved 72 in NLO QCD calculations 

The relevance of resolved 72 within the framework of NLO 
parton level calculations of dijet cross sections in ep col- 
lisions has been investigated using JETVIP |l3|], the only 
NLO parton level MC program including both direct and 
resolved photon contributions. In specifying the powers 
of a and a s corresponding to various Feynman diagrams 
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1 ,4<PX2.4 GeV, 0.05<y<0.95, E T >5 GeV 

- 2.5 




Fig. 2. Upper three plots: diparton 
cross sections, corresponding to target 
■y T and 7^ and using QED as well as 
QCD improved PDF of the latter, plot- 
ted as functions of 77, Et and x~, for 
1.4 < P 2 < 2.4 GeV 2 , 0.05 < y < 0.95, 
Et > 5 GeV, without any restriction on 
r). Lower three plots: ratia of the con- 
tributions of resolved (using PDF of 
|S|) to those of y T (evaluated with PDF 
of integrated over the whole region 
0.05 < y < 0.95, as well as in three in- 
dicated subintervals. 



1.4<P 2 <2.4 GeV 2 , 0.05<y<0.95, E T >5 GeV, -3<7)<0 




Fig. 3. The same as in Fig. |2J, but for 
experimentally motivated restricted re- 



x < gion -3 < 77 < 0. 



we discard one common power of a coming from the ver- 
tex where the virtual photon is emitted by the incoming 
electron. This vertex is also left out in diagrams of Fig. ^|. 

JETVIP contains full set of partonic cross sections for 
the direct photon contributions up the order aa 2 . Exam- 
ples of the corresponding diagrams are shown in Fig. [6|a,b. 
To go one order of a s higher and perform complete calcu- 
lation of the direct photon contributions up to order aa^ 
would require evaluating tree diagrams like that in Fig. 
|e, as well as one-loop corrections to diagrams like in Fig. 
|b and two-loop corrections to diagrams like in Fig. |fj|a. 
So far, such calculations are not available. In addition to 
complete 0(aa 2 ) direct photon contributions JETVIP in- 
cludes also the resolved photon ones with partonic cross 
sections up to the order a^, exemplified by diagrams in 
Fig. ||c,d. The rationale for including in the resolved chan- 
nel terms of the order is discussed in detail in || . Once 



the concept of virtual photon structure is introduced, part 
of the direct photon contribution (which for the virtual 
photon is actually nonsingular) is subtracted and included 
in the definition of PDF of 7*. To avoid misunderstand- 
ing we shall henceforth use the term "direct unsubtracted" 
(DIR uns ) to denote NLO direct photon contributions be- 
fore this subtraction, reserving the term "direct" (DIR) 
for the results after it. In this terminology the complete 
JETVIP calculations are given by the sum of direct and 
resolved parts and denoted DIR+RES. In JETVIP only 
the convolution of QED splitting term (plus some finite 
terms) corresponding to 7^ 

gQ ED (x, P\ M 2 ) = |-3e 2 (x* + (1 - xf) In ^. (5) 

with a 2 s partonic cross sections are subtracted from DIR uns 
calculations. We recall that in any NLO DIR un s calcula- 
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2.4<P 2 <4.4 GeV 2 , 0.05<y<0.95, E T >5 GeV, -3<7)<0 





tion both ^ and 7^ are taken into account exactly up to 
the order aa 2 s . Introducing the concept of resolved 7^ and 
7£ implies the replacement of the convolution (denoted 
er(PSP)) of photon splitting terms (@ for 7J and (H) for 



7£) and order a s partonic cross sections with the contri- 
bution (denoted ut.l (RES)) of the resolved 7^. L . The net 
effect of this operation is thus the addition to cr(DIR uns ) 
of the differences A T , L = °t,l (RES) - cr T , L (PSP) 

a(DIR + RES) = a(DIR uns ) + A T + A L . (6) 

The appropriate measure of the relevance of resolved 7^ 
in NLO calculations is thus the ratio 



r 



NLO 
LT 



(Et,t)) 



A L (E T ,r,) 
A T (E T ,r])' 



(7) 



Note that as for the LO QCD calculations the correspond- 
ing measure is the ratio <tl (RES) /cry (RES), the relevance 



G„(x„r 




b) G,(x,)' 



Fig. 6. Examples of diagrams contributing to dijet production 
in ep collisions at the orders aa s (a), ao? s (b,c), and aal (d,e) 
taking into account that the upper blobs representing quark 
distribution functions of the photon are proportional to a. 



of 7£ in hard collisions is in general different at LO and 
NLO orders. For 7^ the cross section <tl (RES) is given by 
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the convolution of QCD improved PDF of j L with par- 
tonic cross sections up to the order a^. 

To include the effects of resolved j L , we have modified 
JETVIP with the help from Bjorn Potter in three places 
by adding: 

— the flux (|) of 72, 

— the photon splitting term (^) corresponding to 7]?, 

— the call to PDF of 7 £. 

We have checked our modifications against HERWIG as 
well as internally within JETVIP. In the first case we com- 
pared LO JETVIP results for -f L with analogous results 
obtained with HERWIG 5.9 for the same QCD improved 
PDF of initial j L - Taking into account small differences 
between the way JETVIP and HERWIG 

— set the scale of PDF and a s , 

— treat (light) quark mass effects, 

— reconstructs kinematics from generated x 1 , 

we have found very satisfactory agreement in both shape 
and absolute normalization of resulting distributions in all 
three variables 2^,77 and Et- 

The modification of the photon splitting term (||) to 
the form appropriate for 7^ has been checked by compar- 
ing JETVIP results for cr L (PSP) with LO JETVIP results 
in the resolved channel obtained with purely QED expres- 
sion (g) for light quark distribution functions. Apart from 
the opposite sign, the latter should be equal to the former 
as, indeed, it turned out to be the case to within a few %. 

6.1 Hadronization corrections 

Any meaningful comparison of JETVIP results with ex- 
perimental data must involve estimates of the effects de- 
scribing the conversion of partons to hadrons. These hadro- 
nization corrections are not simple to define, but adopting 
the definition used by experimentalists fllif , we have found 
|llj that they depended sensitively and in correlated man- 
ner on the pseudorapidity and transverse energy of jets. 
For E T = 5 GeV, hadronization corrections become large 
and strongly model dependent for 77 < —2.5. We have 
therefore restricted our analysis to —2.5 < rj < 0, where 
they are flat in 77 and do not exceed 10%. 

6.2 Results 

We have redone the calculation of Q using QCD improved 
PDF of 7£, but otherwise with the same assumptions con- 
cerning renormalization and factorization scales ^| and for 
identical kinematical region 

-2.5 < 77 < 0, E^ ] > 7, E^ ] > 5 GeV. 

The resulting distributions da/drj and da/dEx correspond- 
ing to the second window in P 2 are shown in Fig. ^. 
We plot there separately all three contributions on the 

8 In PDF of 7 £ we set A 2 QCU = 0.1 GeV 2 . 
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Fig. 7. Comparison of nontrivial parts At and Al of the con- 
tributions of y T and 7£ to da/drj and da/dEr distributions. 
The results of direct unsubtracted and full calculations using 
in the resolved channel y T only are shown as well. 



r.h.s. of (g), as well as their sum defined in (Q) but in- 
cluding the contributions of 7^, only. Note that both At 
and Al are almost flat in rj and rapidly falling in Et, 
the latter fall-off being faster for At, as expected due to 
harder shape of PDF of j T . The resulting r^°(-Er,ry) 
rises slowly from about 0.2 at 77 = —2.5 to 0.35 at 77 = 0, 
but decreases appreciably with Et- Integrated over Et, 
we find r^Jf ^) ~ 0.3, but for Et close to the lower cut- 
off E T = 7 GeV, this ratio increases to about 0.5. Note 
also that for 77 close to tj ~ 0, Al approaches the results 
of DIR uns calculations. 

Increasing the photon virtuality: 

— reduces the relevance of resolved j T as measured by 
the ratio ZiT/c(DIR uns ), but 
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Inclusion of the contributions of resolved 7^ in phe- 
nomenological analyses of HERA data on dijet produc- 
tion helps bring the theoretical predictions closer to the 
HI data jll), but a thorough analysis of the evidence for 
resolved 7^ in current HERA data requires detailed dis- 
cussion of a number of points, and is beyond the scope of 
this paper. 



7 Summary and conclusions 

We have analyzed the contributions of resolved 7^ to vir- 
tual photon structure function and dijet cross sec- 
tions measured at HERA, using the recently constructed 
parameterization of QCD improved PDF of 7?. The con- 
tributions of resolved 7^ were shown to be nonnegligible 
with respect to those of "f T , but their relevance depends 
on the order of QCD calculations employed and kinemat- 
ical region considered. Within the LO QCD and in the 
kincmatical regions accessible at LEP and HERA, they 
amount typically to 40 — 50% of those coming from re- 
solved j T , but in parts of phase space (small y and x 7 or 
low Et) this number is even larger. Within the NLO cal- 
culations of virtual photon interactions the relative impor- 
tance of 7£ with respect to 7^ is smaller, but still clearly 
of phenomenological relevance. In both cases the effects of 
QCD improved PDF of 7^ are clearly observable. 

We are grateful to J. Cvach, C. Friberg and B. Potter for inter- 
esting discussions concerning the structure and interactions of 
longitudinal virtual photons and to B. Potter for help in mod- 
ifying JETVIP. This work was supported in part by Grant 
Agency of the Academy of Sciences of the Czech Republic un- 
der the grants No. A1010821 and B1010005. 
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Fig. 8. The same as in Fig. | but for 4.4 < P 2 < 10 GeV 2 . 



— increases the relative importance of resolved 7^ with 
respect to resolved j T as measured by the ratio r^jP . 

This is illustrated in Fig. ||, which shows the same plots 
as in Fig. but for 4.4 < P 2 < 10 GeV 2 . In this interval 
the mean value of rf^P is about 0.38, but for Et close to 
E T = 5 GeV it approaches unity. Rising the cut-off E T 
reduces the relevance of 7^ with respect to 7^, for much 
the same reasons as in LO calculations. 
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7^ is determined by two circumstances: the presence of 
"large log" ln(M 2 /P 2 ) in PDF of 7^ and the difference in 
shapes of PDF of j T and 7^. At very large value of the 
ratio M 2 / P 2 the first effect is clearly more important and 
leads to dominance of resolved j T . However, in presently 
accessible range at HERA this "large log" is fairly small 
number around 3 and thus the fact that PDF of j T are 
harder than those of 7^ plays equally important role. 
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